PROCEEDINGS 

0F SCIENCE 



The low-momentum ghost dressing function and the 
gluon mass 



J. Rodnguez-Quintero 

Dpto. Ffsica Aplicada, Fac. Ciencias Experimentales; Universidad de Huelva, 21071 Huelva; 
Spain. 
E-mail: 



jose.rodriguez@dfaie. uhu . es 



We study both regular (the zero-momentum ghost dressing function not diverging), also named 
decoupling, and critical (diverging), also named scaling, Yang-Mills propagators solutions by an- 
alyzing the low-momentum behaviour of the ghost propagator Dyson-Schwinger equation (DSE) 
in Landau gauge, assuming for the truncation a constant ghost-gluon vertex, as it is extensively 
done, and a simple model for a massive gluon propagator. The asymptotic expression obtained for 
the regular or decoupling ghost dressing function up to the order ff{q 2 ) fits pretty well the low- 
momentum ghost propagator obtained through the numerical integration of the coupled gluon and 
ghost DSE in the PT-BFM scheme and, when the size of the coupling renormalized at some scale 
approaches some critical value, the PT-BFM results seems to trend to the the scaling solution as a 
limiting case. 
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1. Introduction 

The low-momentum behaviour of the Yang-Mills propagators derived either from the tower of 
Dyson-Schwinger equations (DSE) or from Lattice simulations in Landau gauge has been a very 
interesting and hot topic for the last few years. It seems by now well established that, if we as- 
sume in the vanishing momentum limit a ghost dressing function behaving as F(q ) ~ (q 2 ) aF and a 
gluon propagator as A(q 2 ) ~ (q 2 ) aa ~ l (or, by following a notation commonly used, a gluon dressing 
function as G(q 2 ) = q 2 A(q 2 ) ~ {q 2 ) a °), two classes of solutions may emerge (see, for instance, the 
discussion of refs. [jl|, ||]) from the DSE: (i) those, dubbed "decoupling", where a,p = and the sup- 
pression of the ghost contribution to the gluon propagator DSE results in a massive gluon propagator 
(see [||, ||] and references therein); and (ii) those, dubbed "scaling" , where (Xf ^ and the low- 
momentum behaviour of both gluon and ghost propagators are related by the coupled system of DSE 
through the condition lap + etc = implying that F 2 {q 2 )G{q 2 ) goes to a non-vanishing constant 
when q 2 — > (see [Q, ^] and references therein). 

Lattice QCD results appear to support only the massive gluon (ag = 1) or scaling solutions 
(see [|7[ ^ 1 , 1C , 1 1, 12] and references therein), and also pinching technique results (see, for instance, 
Jl3| , p| ] and references therein), refined Gribov-Zwanziger formalism (see [|l5|]) or other approaches 



like the infrared mapping of X<p 4 and Yang-Mills theories in ref. fll6| ] or the massive extension of the 
Fadeev-Popov action in ref. [|l7|] appear to point to. 

In the present note, we briefly review the work of refs. ^ ^, |l8| ], where it is established how 
both types of IR solutions for Landau gauge DSE emerge and how the transition between them may 
occur, and that of ref. [19] which extends the previous studies by the analysis of the results [|2C|] 



obtained by solving the coupled system of Landau gauge ghost and gluon propagators DSE within 
the framework of the pinching technique in the background field method [^l|] (PT-BFM) 

2. The two kinds of solutions of the ghost propagator Dyson-Schwinger equation 

As was explained in detail in refs. the low-momentum behavior for the Landau gauge 

ghost dressing function can be inferred from the analysis of the Dyson-Schwinger equation for the 
ghost propagator (GPDSE), which can be written diagrammatically as 
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That analysis is performed on a very general ground: one applies the MOM renormalization pre- 
scription, 



F R {y, 2 ) = M 2 A«(M 2 ) = 1 



(2.2) 



where \i 2 is the subtraction point, chooses for the ghost-gluon vertex, 



ff c (-q,k;q-k) = igof^iqvH^ + iq-k^^k) 



cabc 



(2.3) 
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to apply this MOM prescription in Taylor kinematics {i.e. with a vanishing incoming ghost momen- 
tum) and assumes the non-renormalizable bare ghost-gluon form factor, H\{q,k) =H\,to be constant 
in the low-momentum regime for the incoming ghost. Then, the low momentum-behaviour of the 
ghost dressing function and the gluon propagator is supposed to be well described by 
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and one is finally left with: 
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(2.7) 



It should be understood that the subtraction momentum for all the renormalization quantities is /I 2 . 
The case a F ^ leads to the so-called scaling solution, where the low-momentum behavior of the 
massive gluon propagator forces the ghost dressing function to diverge at low-momentum through 
the scaling condition: 2a F + OCg = (0Cg = 1 is the power exponent when dealing with a massive 
gluon propagator). As this scaling condition is verified, the perturbative strong coupling defined in 
this Taylor scheme [22], a F = g F /(4n), has to reach a constant at zero-momentum, 
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as can be obtained from Eqs.(2.4,2.6). The case a F = corresponds to the so-called decoupling solu- 



tion, where the zero-momentum ghost dressing function reaches a non-zero finite value and eq. ( J2.q ) 
provides us with the first asymptotic corrections to this leading constant. This subleading correction 



is controlled by the zero-momentum value of the coupling defined in eq. (|2.7|), which is an extension 



of the non-perturbative effective charge definition from the gluon propagator [ |23| ] to the Taylor ghost 



gluon coupling [|24[]. As a consequence of the appropriate amputation of a massive gluon propagator, 



where the gluon mass scale is the same Rl-invariant mass scale appearing in eq. (|24|), this Taylor 
effective charge is frozen at low-momentum and gives a non-vanishing zero-momentum value. 



3. Comparison with numerical results from coupled DSE's 



We shall now compare the formulas given by eqs. (2^,2^) with some numerical results for the 
gluon propagator and ghost dressing function obtained by solving the coupled system of gluon and 
ghost DS equations obtained by applying the pinching technique in the background field method 
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(PT-BFM) [21] (see also [|14|] and references therein). In the PT-BFM scheme for the coupled DSE 
system, the ghost propagator DSE is the same of Eq. (2.1), where the approximation H\ = \, and the 
gluon DSE is given by 

(l + Gfo 2 )) 2 ' 



A(q 2 ) 



\S(iv - = q 2 Siiv - q^qv + «X ( a; 



(3.1) 
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where the external gluons are treated, from the point of view of Feynman rules, as background fields 



(these diagrams should be also properly regularized, as explained in [[14]]). The function 1 + G defined 



in ref. [ |25| ] can be, in virtue of the ghost propagator DSE, connected to the ghost propagator Q24[]. 
The coupled system is to be solved, by numerical integration, with the two following boundary 
conditions as the only required inputs: the zero-momentum value of the gluon propagator and that of 
the coupling at a given perturbative momentum, /I = 10 GeV, that will be used as the renormalization 
point. 

Thus, The PT-BFM framework leaves us with an attractive model for gluon and ghost propaga- 
tors providing quantitative description of lattice data [Q, ^6|] and giving well account of their main 
qualitative features: finite gluon propagator and finite ghost dressing function at zero-momentum. 
Futhermore, the coupled DSE system can be solved with different boundary conditions (see be- 
low). In particular, solutions obtained by keeping the zero-momentum value of the gluon propagator 
fixed (see lefthand plots of fig. [j]) while a(jjL 2 = 100 GeV 2 ) is ranging from 0.15 to 0.1817 are avail- 



able [|20|]. these solutions can be confronted to the asymptotical expressions derived in the previous 
section. 

3.1 Decoupling solutions 

Then, as the gluon propagator solutions in the PT-BFM scheme result to behave as massive ones, 



the eqs. ( |2.4| , [2.6[ ) must account for the low-momentum behaviour of both gluon propagator and ghost 
dressing function with Hi = 1 and Ofr(0) given by eq. (2.7), with 0lr(/l 2 ) = g R (/x 2 ) / (4n) being 
fixed, as a boundary condition for the numerical integration of the coupled DSE for each particular 
solution of the family (see tab. [j]). Furthermore, the zero-momentum values of the ghost dressing 
function, Fr(0) and of the gluon propagator, Ar(0), can be taken from the numerical solutions of 
the DSE (for any value of the a{pL = 10GeV)). These altoghether with the gluon masses obtained 
by the fit of eq. ( 2.4 ) to the numerical DSE gluon propatator solutions (see the left plot in fig. [j], for 
a(pi) = 0.16, and the results for a(ju) = 0.15,0.16,0.17 in tab. [IJ, taken from ref. Jl^]), provide us 
with all the ingredients to evaluate, with no unknown parameter, eq. ( |2.6[ ). 

Indeed, the expression given by eq. ( ^^ ) can be succesfully applied to describe the solutions all 
over the range of coupling values, cc(jl), at pL = 10 GeV (provided that they are not very close of 
the critical coupling that will be defined in the next subsection). This can be seen, for instance, for 
a = 0.16, in the right plots of fig. [I] and it is shown for a = 0.15,0.16,0.17 en ref. [Oj. 
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A(q>=10) 



oe(10 GeV) = 0.16 



F(q",H=10) 



a(10GeV)=0.16 




q" (GeV ) 



q" (GeV ) 



Figure 1: Gluon propagator (left) and ghost dressing function (right) after the numerical integration of the 
coupled DSE system for a(fJ, = lOGeV) = 0.16 taken from [ |20| ] . The curves for the best fits to gluon 
propagator and ghost propagator data explained in the text appear as red dotted lines, the same for the black 
dotted line in the lefthand plot but retaining only the logarithmic leading term for the asymptotic ghost dressing 
function by dropping the —11/6 away. In the righthand plots, the red dotted lines correspond to apply. 





«r(0) 


M (GeV) [gluon] 


0.15 


0.24 


0.37 


0.16 


0.30 


0.39 


0.17 


0.41 


0.43 



Table 1: Gluon masses and the zero-momentum non-perturbative effective charges, taken from ref. 
obtained as discussed in the text. 



and 



3.2 The "critical" limit 

There appears to be a critical value of the coupling, a cr it = oc(/J. 2 ) ~ 0.182 with pt = 10 Gev, 
above which the coupled DSE system does not converge any longer to a solution [|20|]. As a matter 
of the fact, we know from refs.j^, [19|] that the scaling solution implies for the coupling 



Or, 



47T 2N c A 2 (n 2 )B{n 2 



(3.3) 



where B(/i 2 ) is determined by the gluon propagator solution that is supposed to behave as eq. (2.4), 
and A(/i 2 ) by the ghost propagator that should behave as eq 



L6[ ) in the case ap / 0. Again, /i 2 is the 
momentum at the subtraction point. This is also shown in ref. [[I]], where only the ghost propagator 
DSE with the kernel for the gluon loop integral is obtained from gluon propagator lattice data. In 
the analysis of ref. [|J], a ghost dressing function solution diverging at vanishing momentum appears 
to exist and verifies eqs. (2^,3^), while regular or decoupling solutions exist for any a < a cr i t . In 



ref. [|19|], a more complete analysis is performed by studying again the dressing function computed 

a(n 2 ), at jx 2 = 100 GeV 2 @. 



by solving eq. ( |3.1| ) for the different values of the coupling, a 
The ghost dressing function at vanishing momentum, F(0,/i 2 ), is described by the following power 
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behaviour, 

F(0) ~ (a crit -a(At 2 ))-^ 2) , (3.4) 

where k:(jU ) is a critical exponent (depending presummably on the renormalization point, /I 2 ), sup- 
posed to be positive and to govern the transition from decoupling (a < Gkrit) to the scaling (a = Gfcrit) 
solutions; and where we let a cr jt be a free parameter to be fitted by requiring the best linear correla- 
tion for log[F(0)] in terms of log[a cr i t — Of]. In doing so, the best correlation coefficient is 0.9997 for 
CUcrit = 0.1822, which is pretty close to the critical value of the coupling above which the coupled DSE 
system does not converge any more, and K"(/x 2 ) = 0.0854(6). This can be seen in fig. |2|(a), where 
the log-log plot of Fg(0) in terms of otrit — en is shown and the linear behaviour with negative slope 
corresponding to the best correlation coefficient strikingly indicates a zero-momentum ghost prop- 
agator diverging as a — > a c rit- Nevertheless, no critical or scaling solution appears for the coupled 
DSE system in the PT-BFM, although the decoupling solutions obtained for any a < a cl -i t = 0.1822 
seem to approach the behaviour of a scaling one when a — > a cr i t . The absence of the scaling solution 



in the PT-BFM approach can be well understood by analysing eq. (3. 1 ) as explained in ref. [ 19] 



A(qV) 




Figure 2: (a) Log-log plot of the zero-momentum values of the ghost dressing function, obtained by the 
numerical integration of the coupled DSE system in the PT-BFM scheme, in terms of a cr it — 05. a = a(jtt = 
lOGeV), the value of the coupling at the renormalization momentum, is an initial condition for the integration; 
while a cr i t is fixed to be 0.1822, as explained in the text, by requiring the best linear correlation, (b) Gluon 
propagator solutions in terms of q 2 for the same coupled DSE system for different values of a(n = lOGeV), 
all very close to the critical value, ranging from 0.18 to 0.1817 . 

When approaching the critical value of the coupling, the gluon propagators obtained from the 
coupled DSE system in PT-BFM must be also thought to obey the same critical behaviour pattern as 
the ghost propagator. In the PT-BFM, the value at zero-momentum being fixed by construction M, 



20], one should expect that, instead of decreasing, the gluon propagator obtained for couplings near 



to the critical value increases for low momenta: the more one approaches the critical coupling the 
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more it has to increase. This is indeed the case, as can be seen in fig. ^b). This also implies that, 
near the critical value, the low momentum propagator does not obey eq. (^4j) and that consequently 
eq. (2.6) does not work any longer to describe the low momentum ghost propagator. 



4. Conclusions 

The ghost propagator DSE, with the only assumption of taking H\ (q, k) from the ghost-gluon 



vertex in eq. (2.3) to be constant in the infrared domain of q, can be exploited to look into the low- 
momentum behaviour of the ghost propagator. The two classes of solutions named "decoupling" and 
"scaling" can be indentifled and shown to depend on whether the ghost dressing function achieves 
a finite non-zero constant (ap = 0) at vanishing momentum or not {dp ^ 0). The solutions appear 
to be dialed by the size of the coupling at the renormalization momentum which plays the role of a 
boundary condition for the DSE integration. 

We applied a model with a massive gluon propagator to obtain the low-momentum behaviour 
of the ghost propagator that results to be regulated by the gluon mass and by a regularization- 
independent dimensionless quantity that appears to be the effective charge defined from the Taylor- 
scheme ghost-gluon vertex at zero momentum. Then, we demonstrated that the asymptotic de- 
coupling formula (ap = 0) successfully describes the low-momentum ghost propagator computed 
trhough the numerical integration of the coupled gluon and ghost DSE in the PT-BFM scheme. We 
also show that the zero-momentum ghost dressing function tends to diverge when the value of the 
coupling dialing the solutions approaches some critical value. Such a divergent behaviour at the 
critical coupling corresponds to a scaling solution where, if the gluon is massive, ar = —1/2. 
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